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Abstract. We examine the nucleoli's electromagnetic form factors in a Poincare-covariant Faddeev framework. The three- 
quark core contributions to the form factors are obtained by employing a quark-diquark approximation. We implement the self- 
consistent solution for the quark-photon vertex from its inhomogeneous Bethe-Salpeter equation. We find that the resulting 
transverse parts which add to the Ball-Chiu vertex have no significant impact on nucleon magnetic moments. The current- 
quark mass evolution of the form factors agrees with results from lattice QCD. 
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Investigating the nucleon's electromagnetic structure continues to be an experimental challenge in contemporary parti- 
cle physics. High-precision electron scattering measurements have provided information on the nucleon's electromag- 
netic form factors as well as longitudinal vs. transverse momentum fractions and the distribution of spin and orbital 
angular momentum among its constituents; see [1, 2] for reviews. Approaching these issues from a theoretical per- 
spective aims at an understanding of the hadrons' substructure in terms quarks and gluons, the fundamental degrees of 
freedom in Quantum Chromodynamics. 

In a Dyson-Schwinger/bound-state equation approach (see [3-5] for reviews) baryons are obtained as solutions of 
the covariant Faddeev equation. It describes the binding of a baryon by iterated quark-quark correlations and thereby 
constitutes the three-body analogue of a quark-antiquark Bethe-Salpeter equation (BSE). The nucleon's Faddeev 
equation was recently solved in a numerical setup where the full covariant structure of the nucleon amplitude was 
implemented [6]. The quark-quark kernel was modeled by a rainbow-ladder interaction, i.e. a dressed gluon exchange 
between any two quarks, where the quark propagator was consistently determined from its Dyson-Schwinger equation 
(DSE). A rainbow-ladder truncation excludes the presence of pion-cloud corrections which have significant impact on 
the chiral structure of hadrons. It therefore reflects the properties of a hadronic 'quark core'. 

The computation of electromagnetic form factors in this setup requires knowledge of another ingredient, namely the 
dressed quark-photon vertex. A frequently used truncation of the Faddeev equation is the quark-diquark model, where 
the binding of baryons is traced back to the interaction between quarks and effective scalar and axial- vector diquarks 
within a baryon [7, 8]. It was demonstrated in Ref. [6] that this simplification only introduces a <~ 5% change in the 
nucleon mass. Hence, as an intermediate step towards a fully self-consistent determination of electromagnetic form 
factors from the covariant Faddeev equation, we investigate the impact of the structure of the quark-photon vertex on 
nucleon form factors in the quark-diquark approach. 



The nucleon's electromagnetic current is expressed in terms of two form factors: the Dirac and Pauli form factors 
Fi(Q 2 ) and F2(Q 2 ), or the Sachs form factors as their linear combinations Ge = F\ — xFi and Gm = F\ + F2, with 
T = Q 2 /{AMji). Their static values represent the proton and neutron charges X p - n = (1,0), (anomalous) magnetic 
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INTRODUCTION 



NUCLEON ELECTROMAGNETIC FORM FACTORS 




6F/(0), 
6G' E (0) 



Kr\= -6F 2 '(0), 
Md=-6GV(0). 



(1) 



Electromagnetic current conservation provides a construction principle for the current in a given model frame- 
work [9]. In the quark-diquark approach this construction involves quark and diquark impulse-approximation di- 
agrams, but also a coupling to the exchanged quark between quark and diquark as well as seagull terms, i.e. the 
photon's coupling to the diquark amplitudes. Upon resolving the diquark-photon vertex into its constituents, all pho- 
ton couplings are specified by two quantities: the quark-photon vertex and the seagull vertices. The details of this 
construction are described in Refs. [8, 10] for the nucleon and Ref. [11] for the A-baryon, respectively. 

Herein we want to investigate the structure of the quark-photon vertex in the context of electromagnetic form factors. 
The most general form of the vertex can be expressed by a sum of the Ball-Chiu term [12] and a purely transverse 
contribution: 

it Q) = [»V£a + 2 ^ ('* a a + As) ] + r£ (*, Q) • (2) 
The Ball-Chiu part involves the quantities 
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where k± — k±Q/2 are the in- and outgoing quark momenta and A(k 2 ) and B(k 2 ) = M{k 2 )A{k 2 ) are the dressing 
functions of the quark propagator S(k) = A~ l (-it + M) / (k 2 +M 2 ). Eq.(2) satisfies the Ward-Takahashi identity 
Q^Fq(k, Q) = S^ 1 (k + ) — S~ l (k-). The transverse part can be written as 

-zlt = fi $ + ifi Y$<2 + h h k Q [# , *] + U tf*T <S + % (if s +f 6 k-Q® + fik + if & k T @), (4) 

where the fi(k 2 , k ■ Q, Q 2 ) are scalar dressing functions. We used the abbreviations jj = Tq" ' y v ', kj = Tg V k v and 

yf = T^jj, with Tq v = d^ v — Q^Q V being a transverse projector with respect to the photon momentum Q. 
The quark-photon vertex is obtained from its inhomogeneous Bethe-Salpeter equation [13], given by 

[r${k,Q)-z 2 i<f] afi - \ f p K aa ,p, p [s(k' + )r^(k\Q)s(k'_)] a , pi , (5) 

where Z2 is the quark renormalization constant and K aa ipip the rainbow-ladder kernel that also appears in the quark 
DSE and the homogeneous meson and diquark BSEs. It reads 

W=^ 4 -^«.^, (6) 

with k being the gluon momentum, and involves an effective interaction a(k 2 ) which defines the model input of the 
approach. We use the parametrization of Ref. [14] for a(k 2 ) whose infrared term is modeled by two parameters: an 
infrared scale A and a width parameter 77 . Adjusting the scale A to reproduce the experimental pion decay constant 
yields a good description of pseudoscalar-meson, vector-meson, nucleon and A ground states. Inflating this scale 
towards the chiral limit mimicks the properties of an effective quark core: it generates hadron masses which are 
overestimated in the chiral region and must be dressed, for instance, by pionic corrections; see Refs. [11, 15] for 
a further discussion. These studies exhibited only a modest dependence on the width parameter 77; i.e., the quoted 
observables are not sensitive to the details of the interaction in the infrared. 

To solve the inhomogeneous BSE it is convenient to express the vertex in an orthonormal basis, for instance 
defined by the elements rf (k, Q) G {y? A/2, k T ^ , g M } x {1 , @ , kj , It 0} which satisfy the orthogonality relation 
Tr{tf Tj'} = 4dij. Since the inhomogeneous equations for the transverse and longitudinal parts decouple, it is sufficient 
to consider the eight transverse elements alone. The longitudinal result reproduces the longitudinal projection of 
the Ball-Chiu vertex and thereby the Ward-Takahashi identity, and purely longitudinal terms do not contribute to 
nucleon form factors because of current conservation: Q^JP = 0. The inhomogeneous BSE self-consistently generates 
a timelike vector-meson pole in the quark-photon vertex at Q 2 = —m 2 p which significantly increases the charge radii 
of pseudoscalar and vector mesons [13, 16]. 



RESULTS AND DISCUSSION 



Upon implementing the numerical solution for the quark-photon vertex according to Eqs. (5-6) in the nucleon form 
factor setup, one can quantify the impact of the purely transverse structure (4) of the vertex. In previous form factor 



TABLE 1. Static electromagnetic properties of the nucleon at the physical pion 
mass m K = 140 MeV. The magnetic moments jl p ' n are dimensionless and the 
squared charge radii R E '" M = (t^ n M M^) 2 are given in (GeVfm) 2 . The first row 
shows the results obtained with a Ball-Chiu vertex augmented by the ansatz (7); 
the second row those with the quark-photon vertex solution from the inhomoge- 
neous BSE. The brackets indicate the sensitivity to the model parameter rj. The 
last row collects experimental results [17, 18]. 





I M p 




R p E 




Km 




BC + p 
Full vertex 


2.56(6) 
2.52(12) 


-1.59(3) 
-1.59(7) 


0.55(2) 
0.58(1) 


0.00(1) 
0.01(1) 


0.46(1) 
0.50(1) 


0.45(1) 
0.49(1) 


Exp. 


2.79 


-1.91 


0.69(2) 


-0.10 


0.64(5) 


0.67(2) 



studies this transverse part was modeled by an effective p— meson contribution [10, 19]: 

r$ = -^-^e-^ 2 H)r vc (k,Q), (7) 

gp Q 2 + mj 

where T vc (k, Q) is the calculated vector-meson amplitude (see Ref. [19] for details). In consistency with the behavior 
of the pion's charge radius, the Ball-Chiu part and Eq. (7) each roughly contribute ~ 50% to the nucleon's squared 
charge radii {r p E ) 2 and (f^") 2 , whereas Eq. (7) cancels the Ball-Chiu piece in the neutron's electric charge radius (r E ) 2 . 
On the other hand, this ansatz does not contribute to the nucleon magnetic moments due to the inherent suppression 
with Q 2 . The resulting static electromagnetic properties allow for the interpretation as a quark core which agrees well 
with lattice data at larger quark masses where meson-cloud effects are no longer important. 

It is visible from Table 1 that the full vertex-BSE solution does not significantly change this behavior. The resulting 
nucleon charge radii differ by < 5% from those obtained via Eq. (7). The fractions of jj, p and jj," which come from 
the transverse vertex-BSE solution of Eq. (4) owe to the amplitudes fi and f\. Those two substantially contribute to a 
quark's anomalous magnetic moment [20] where they roughly cancel the Ball-Chiu part. Nonetheless, they have only a 
small impact upon the nucleon magnetic moments which are still sufficiently well represented by the Ball-Chiu ansatz 
alone. Moreover, we have verified that the Q 2 -evolution of all electromagnetic form factors is essentially identical for 
both choices of transverse contributions up to Q 2 <~ 3 GeV 2 . 

In Fig. 1 we show the current-mass evolution of the nucleon magnetic moments and charge radii. The left panels 
depict the dimensionless magnetic moments jj, pn and squared charge radii R E " M := {^"mMn) 2 which exhibit a 
moderately rising behavior with the pion mass. In the remaining panels we show the isovector and isoscalar anomalous 
magnetic moments and isovector Dirac and Pauli radii as obtained from the form factor combinations F>' = Ff — F" 
and Ff = Ff +F", with i = 1,2. To allow for a sensible comparison of dimensionful charge radii and magnetic 
moments in static nuclear magnetons with lattice results we plot the following rescaled quantities: 

^ = K %> ^, = r 2 {^)\ with MfKf=M? + (^)V/k)), (8) 



where K, r and Mn are our calculated results and M* ef , with f{m\ ) = 0.77/(1 + (mjO.65 GeV) 4 ) and M Q = 0.9 GeV, 
is a reference mass which describes the lattice results for Mn of Refs. [21-26] reasonably well. 

Fig. 1 shows a good agreement between our calculated magnetic moments and the lattice data. This is especially 
relevant at higher pion masses where the pion-cloud dressing effects of the nucleon are diminished. In combination 
with the observed missing chiral curvature, our results can be viewed to represent the quark core of the nucleon. The 
implementation of such chiral effects, together with a form factor determination beyond the quark-diquark model, 
remains a task for future investigations. 
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FIGURE 1. Static electromagnetic properties of the nucleon as functions of the squared pion mass. Left panels: dimensionless 
magnetic moments and squared charge radii in units of (GeVfm) 2 . Middle and right panels: isovector and isoscalar anomalous 
magnetic moments and squared isovector charge radii, rescaled as described in the text and compared to a selection of recent lattice 
results [21-26]. The bands show the dependence on the model parameter tj; stars denote experimental values [17, 18]. 
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